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EFFECT STRUCTURE AND FOUNDATION INTERACTION 


Daniel ASCE, and Nathan Newmark,2 ASCE 


SYNOPSIS 


the design structures subjected shock-type loadings, such earth- 
quakes ground shock from nuclear blasts, the shock described 
record the foundation motion generally obtained from structure that has 
been subjected similar type shock. This often presented 
the form response spectra. When this information used design 
analyze structure, the assumption made that the spectrum obtained from 
the foundation motion particular structure accurate representation 
the effects the shock. This implicitly assumes that the dynamic response 
the structure does not affect the recorded foundation motion. 

Analytical models for particular types structures are developed and 
analyzed this paper toward the purpose studying the type and the magni- 
tude the errors that may occur when the record the foundation motion 
specific structure used determine the response another dissimilar 
structure. The effect various assumed structural and foundational para- 
meters these errors examined. 


INTRODUCTION 


many engineering problems necessary estimate the required 


Note.—Discussion open until March 1962. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the Engineering Mechanics Division, Proceedings the American Society Civil 
Engineers, Vol. 87, No. October, 1961. 

Capt., Corps Engrs., U.S. Army Engr. Waterways Experiment Sta., Vicksburg, 
Miss. 

Head, Dept. Civ. Engrg., Univ. Urbana, 
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shocks. Such problems arise shock-mounting sensitive equipment, the 
resistance structure earthquake shock and, more recently, the re- 
sistance structure ground shock from atomic explosion. Also, the 
problem has become important the shock isolation missile systems and 
underground protective structures. The specification the type and intensi- 
ties shocks structure must resist great importance these prob- 
lems. the specification description the shock motion inaccurate, 
the structure may overdesigned underdesigned. Such error may re- 
sult the failure the structure, mayresult structurethat made 
many times stronger and more costly than necessary. 

the design structures subjected shock, the shock generally de- 
scribed record motion another structure that has been 
subjected similar shock. Thus, the assumption the foundation 
motion particular structure accurate representation the shock. 
This implies that the motion the structure does not affect the foundation 
motion, which incorrect; however, this effect has usually been considered 
negligible. The object this paper study the errors that may occur when 
the record the foundation motion particular structure used deter- 
mine the response other dissimilar structures, and analyze the effect 
several structural and foundational parameters these errors. 

not intended use specific types shocks specific structures 
this study, but rather use broad approach order obtain information 
that will lead general understanding the errors that may result. Be- 
cause most the work inthis problem been concerned with aseismic 
design, the examination presented herein will refer building founded 
ground. However, the generalizations used permit this study apply other 
types vibrating systems other types foundations. 

Notation.—The letter symbols adopted for use this paper are defined 
where they first appear, the illustrations the text, and are arranged 
alphabetically, for convenience reference, Appendix 


STRUCTURAL RESPONSE 


Development Problem.—In the determination the response struc- 
tures dynamic load convenient use the so-called response spec- 
strum concept. This concept was described connection with engineering 
seismology Biot 1943 (3), and extensive work has been done since 
may defined plot the maximum value displacement, velocity, 
acceleration single-degree-of-freedom system function the fre- 
quency the system. These maximum values are obtained for given type 
input. the study the response buildings earthquakes, the spec- 
trum obtained, briefly, follows: record obtained during earth- 
quake the motion the base building means accelerometer 
displacement meter installed the basement the building. Using shaking 
table or, more usually, electrical analog computer digital computer, 
this motion record used input the base motion for large number 
fixed-base, single-degree-of-freedom, spring-mass systems, each witha 


Numerals parentheses refer corresponding items the Bibliography, Appen- 
dix 
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different frequency. The maximum response each spring-mass system 
obtained and plotted function the frequency. Shock spectra are ob- 
tained this manner for large number recorded earthquakes and 
envelope these spectra used design spectrum. 

For multi-degree-of-freedom system the response for each mode can 
obtained from the spectrum, and combining the modal responses, the total 
response the system can approximated. The response obtained as- 
suming that the maximum response each mode occurs the same time 
gives upper bound the actual maximum response. 

The development generalresponse spectrum motion recorded 
the basement building implies the assumption that the recorded motion 
accurate description the ground motion. other words, the assump- 
tion made that the vibrating structure does not affect the input foundation 
motion. The question this effect structure-foundation interaction has 
been raised several investigators (4),(8),(13),(20),(47), but has only been 
since 1958 that studies this effect have been made 

There evidence that analysis using the response spectrum concept 
may give results that are overconservative. Tung and Newmark (46), 
using the same method analysis was used for the fixed-base systems 
this study, computed values dynamic shear ten-story building that were 
considerably greater than the design values. Calculated values considerably 
higher than design values have also been reported other authors (1),(4),(11), 
(33),(45). There are many factors that may contribute the great variation 
between computed shears andactualobserved shears. excellent discussion 
some these factors given John Blume, ASCE (4). 

Generally, the discrepancies between observed and computed values re- 
sponse are attributed inadequate recognition the structural properties 
building materials the manner which buildings are fabricated (4),(11), 
however, O’Hara has presented evidence (30),(31) that part these 
discrepancies may related the method which the shock intensity 
described the shock spectrum. The results experimental work (30) indi- 
cate that there error the envelope the individual spectra 
determine the design spectrum. The analyses performed inthis study indicate 
that even using individual spectrum may lead errors. 

General Assumptions.—The response multistory building dynamic 
loading complex because the distribution mass, stiffness each story, 
height each story, and type and quality construction will determine the 
vibration characteristics the structure. The structure has infinite num- 
ber degrees freedom and complete solution tothe problem not prac- 
tical. However, making certain assumptions, complicated structure can 
simplified equivalent dynamic system with afinite number degrees 
freedom. The simplification structure for the purpose making dy- 
namic analysis requires great deal judgment order that the results 
the analysis will closely approximate the response the actual structure. 

flexible structure, such rigid frame building, lateral forces are 
transmitted mainly through bending the columns, whereas stiff struc- 
ture, such shear wall structure, lateral forces are transferred through 
both bending and shearing action the shear walls. For many types struc- 
tures, usually assumed that the building will deflect shear-beam, 
that is, the floors are assumed move parallel each other with the relative 
deflection any story being dependent only the shear that story. These 
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structures are called shear-beam buildings. The mass each floor anda 
portion the mass the columns attached each floor are assumed 
concentrated each floor level. Thus, structure with infinite number 
degrees freedom reduced system which has many degrees 
freedom has concentrated masses. The shear-beam building shown 
Fig. 1(a) can, therefore, represented spring-mass system shown 
Figs. 1(b) and 1(c). The masses represent the lumped mass each floor 
and its associated columns, the springs represent the resistance each story 
shear forces, and the dashpots represent damping effect caused fric- 
tional forces the vibrating building. four-story, shear-beam building was 
used most the analytical studies described this paper. The column 
shears were assumed proportional the relative horizontal displace- 
ments between stories and the spring constants, were assumed remain 
constant regardless the amount deflection, that is, the structure was as- 
sumed deflect linearly elastic manner. Damping was viscous type, 
and was taken either proportional the absolute velocity each story 
the relative velocity between adjacent stories. 

Governing Differential study two different types foun- 
dation conditions are assumed. four-degree-of-freedom, spring- 
mass system assumed attached movable “free” base mass 
which subjected applied, time-dependent force, P(t). This system, 
without damping, shown Fig. 2(a). set simultaneous, second-order, 
differential equations then obtained writing the equation motion for 
each the masses. The general form the equations, which depends the 
type damping considered, is, 


and 
+ + kj (x; Xj-1) = (Xi+1 x; ) =0 (i = 2, 3, 4). .(1b) 


for absolute damping and 


and 


for relative damping. 

the other case, the spring-mass system assumed attached 
rigid “fixed” base and the base undergoes time-dependent displacement, 
This system, undamped, shown Fig. 2(b). The equation motion 
for each mass obtained before and the resulting set simultaneous 
differential equations the same except for the deletion the equation for 
the base mass. both cases only displacements and forcés parallel the 
long axes the spring-mass systems are considered. This means that the 
actual structure only the horizontal component the displacements are being 
considered. 

Numerical Method Solution.—An exact solution the simultaneous dif- 
ferential equations for multi-degree-of-freedom system usually imprac- 
tical, and numerical methods are used. exact solution the equations 
implies the determination the displacement and velocity function 
time. numerical solution consists obtaining numerical values dis- 
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placement and velocity discrete times. addition, numerical-integration 
methods are particularly well suited digital-computer operations. There 
are many different procedures available for the numerical integration the 
differential equations. The method used this study for obtaining the re- 
sponse the spring-mass systems successive approximation method 
integration (24),(27) with assumed linear variation acceleration each 
time interval value 1/6). The time intervals used the numerical 
calculations were chosen that the value was within the stability 
and convergence limits (27). Calculations were made the ILLIAC the 
University Urbana, and the IBM 650 computer the Water- 
ways Experiment Station, Vicksburg, Miss. 

Natural Frequencies.—One the principal parameters defining the dy- 
namic response multi-degree-of-freedom system the set natural 
frequencies. When system freely moves with simple harmonic motion, this 
motion called natural normal mode vibration, and the frequency 
this motion the natural frequency. For multi-degree-of-freedom system, 
the number modes and frequencies equal the number degrees 
freedom the system. For moderately damped systems, the damped and un- 
damped natural frequencies differ only slightly, and customary com- 
pute the undamped natural modes and frequencies. 

For the four-degree-of-freedom systems studied, the equations motion 
for each mass were written the same manner before except that for free 
vibration there external disturbing function. The eigenvalues and the 
corresponding natural frequencies were determined setting the determinant 
the coefficients this set equations equal zero. 

For free system the lowest frequency corresponds rigid body motion 
the system; that is, degenerates zero. Therefore, although there are 
five natural frequencies for the five-mass free systems, only the four non- 
zero values are interest and the lowest nonzero value referred the 
first natural frequency. 


METHOD ANALYSIS 


objective this study todetermine the shock spectra 
computed from the recorded base motion structure can used predict 
accurately the response another dissimilar structure. four-story shear 
building used order introduce higher modes, yet not require extensive 
time for computing the response. assumed that certain mass soil 
with the base the building when the structure subjected ground 
shock. Therefore, two spring-mass models are used represent the struc- 
ture. One “free-base” system representing the structure and the sur- 
rounding soil, and the other “fixed-base” system representing the type 
system generally used dynamic analyses. 

The procedure used compare the accuracy the spectrum method 
predicting response shown diagrammatically Fig. and follows: 


dynamic load applied the base mass series free-base, 
systems. The parameters the spring-mass systems are 
chosen that range frequencies obtained. 

The response each free-base system the dynamic load com- 
puted. 
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record the base displacement, function time, the base 
mass one the free-base, spring-mass systems subjected the dynamic 
load obtained. 

The record the base displacement obtained step used 
input base displacement for series fixed-base, spring-mass systems. 
The fixed-base systems have the same properties the free-base systems 
used step with the exception the base mass. 

The response each fixed-base system the base displacement 
computed. 

The responses computed and are plotted function the fre- 
quency the systems and the results are compared. 


The primary purpose of.a dynamic analysis determine the maximum 
stresses the structure. is, therefore, important determine the maxi- 
mum relative displacement between stories because the stresses are propor- 
tional this quantity. Thus, the term “response” used hereafter will indi- 
cate the absolute value the maximum relative displacement, introducing 
the notation 


Spring-Mass illustrates the various spring-mass systems 
studied. Each system identified two-part code designation. The alpha- 
betic character the designation indicates specific value for each mass 
the spring-mass system and particular distribution spring stiffnesses. 
The numeric characters indicate specific values for the spring stiffnesses. 
The spring stiffnesses are varied for each spring-mass system within par- 
ticular alphabetic group order obtain range frequencies. The code 
tain system. Thus, indicates the absolute displacement the ith mass 
system Al. 

The mass and spring stiffness ratios are maintained constant during par- 
ticular study that the mode shapes each system are the same. the 
and systems the story mass and spring stiffness ratios are all equal 
unity, whereas the systems, the upper two stories have value story 
mass and spring stiffness equal one-half those the lower two stories. 
small number calculations were made for the systems that had two 
equal masses. 

Three different values are assumed for the magnitude the base mass. 
equal the sum the story masses the and systems; one 
the story masses inthe totentimes one the story masses 
the systems. The magnitude the mass soil that acts with the foun- 
dation structure will depend variety factors such the type 
foundation, that pile, footing raft; the geologic characteristics the 
site; the presence water; the physical properties the soil; the amount 
restraint the soil mass; the type failure planes within the soil; and the 
type motion that the structure undergoes. impossible define ac- 
curately this mass soil (17),(19),(42),(44) and the values chosen herein are 
used merely obtain variety conditions. 

Another problem determining suitable model for the soil mass acting 
with the foundation building the character the resistance the soil. 
considering foundation motion, investigators have used many methods 
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define soil stiffness. Various types stiffness constants with large range 
values have been presented (3),(10),(17),(35),(38). Also studies have been 
made propagation stress waves soils and various modeling schemes 
have been used. Smith, A.M. ASCE, (37) used Maxwell model, Voigt 
model and modified Voigt model, and alsodeveloped more complex model. 
Schiffman, ASCE, (36) describes the use Bingham models and 
Burgers models. addition the complexity these types models, the 
problem modeling complicated the lack knowledge about the be- 
havior soils under dynamic loads. Because the advent nuclear wea- 
pons and the associated requirement for underground protective structures 
resist blast-type loadings, there has been increase research deter- 
mining the dynamic properties soils (19),(36),(41),(49). the present time 
(1961) the exact properties are pertinent its dynamic response 
are uncertain. 

Because the many uncertainties soil properties and the lack infor- 
mation about the changes inthese properties under dynamic loadings, and be- 


SYSTEMS 


cause not the purpose this study investigate these effects per se, 
the resistance the soil deformation represented single dashpot. 
Kawasumi and Kanai indicate (18) that studies actual buildings show 
that the vibration energy building dissipates into the ground manner 
“exactly similar that usual viscous damping.” 

Although the problem nonelastic deformation structures important 
(11),(40), not considered this study. The structure assumed re- 
main linearly elastic, which considerably simplifies the numerical calcula- 
tions. However, the application the integration method used can applied 
nonlinear behavior (27). 

The assumed values the masses each alphabetic spring-mass system 
and the spring stiffness distribution are shown Fig. The values the 
fixed-base natural frequencies, and the free-base natural frequencies, 
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were computed for all the systems. The subscripts the and notations 
refer the first (that is, lowest) frequency, second frequency, and on. 

Damping Conditions.—Damping buildings may due partly external 
friction between structural members connections, partly internal fric- 
tion the material, and partly damping from air other adjacent build- 
ings. not known just what type damping most representative 
structural damping (28). For mathematical simplicity, usually assumed 
that the damping viscous damping which taken either proportional 
the absolute velocity each mass, the relative velocity each mass 
with respect adjacent masses, both. Stevens (39) uses model 
based interfloor dashpots; Tung and Newmark (46) use absolute viscous 
damping, whereas Merritt and Housner (22) recommend combi- 
nation both types. this study both types were considered separately. 

customary define the damping system (c) proportion (r) 
the coefficient critical damping For linear, spring-mass-dashpot, 
single-degree-of-freedom system, with circular natural frequency 


For ease identification herein, the notation used designate the pro- 
portion damping applied the base mass, Mo, which proportional the 
absolute velocity used designate the proportion damping ap- 
plied the ith mass, which proportional the absolute velocity 
and used designate the proportion damping which proportional 
the relative velocity between masses and Damping the base mass 
considered proportional the circular natural frequency the first free- 
base mode, and damping the story masses considered proportional 
the circular natural frequency the first fixed-base mode, Thus, 


and 
for systems with absolute damping, and 
and 


for systems with relative damping. 

The amount damping that can expectedin structures not easily de- 
termined. Merit White, ASCE, (48) result tests reinforced- 
concrete storage building, indicates value absolute damping approxi- 
mately critical. Results tests four-story concrete building also 
show (14) damping the order 0.07. Values damping other real build- 
ings have been found (22) range from 0.07 0.40 critical damping. 
this study values damping 10% and 20% were used for the base mass, and 
values and 10% were used for the story masses. higher value damp- 
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ing was used for the base mass considering the fact that the entire dissipation 
energy into the soil represented single dashpot. Calculations re- 
sponse were made using all possible combinations the values damping. 
The values damping previously indicated were used for both absolute and 
relative damping. 

Applied Force Functions.—Four time-dependent force functions were arbi- 
trarily selected for the force applied the base mass. Two these force 
functions, and are shown diagrammatically Fig. The first force 
function, was single-cycle sine pulse with period equal terms 
the time parameter, used. The second force function, was single- 
cycle sine pulse with period These two forces were selected order 
compare the effect the period onthe response the systems. 


(t) 
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FIG, 4.-APPLIED FORCE FUNCTIONS. AND 


After the initial studies with the simple sine pulses, random force func- 
tions, and were developed. Force was obtained making the force 
and time scales proportional the accelerogram the 82° component 
the Vernon, California, earthquake October 1933. The portion the ac- 
celerogram that was used represents actual time sec. Because, 
this study, this value time made equal 100 terms the time para- 
meter unit value equals 0.05 second. this time proportionality 
applied the spring-mass systems their fixed-base periods range 
from 0.13 sec 3.0 sec. Therefore, this particular time scale makes the 
system periods compatible with the range periods for which 
spectra have been obtained (approximately 0.2 sec 3.0 seconds) (9). 

Force selected arbitrarily represent forcing function nearly 
random possible, based studies number earthquake accelero- 
grams. 
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Dimensions.—The maximum value each applied force function ex- 
pressed terms dimensionless quantity multiplied force, Other 
quantities are expressed nondimensional numerical quantities multiplied 
factors related and my. The dimensional relations for the various 
multiplying factors used are follows: 


acceleration, 
displacement, (L) 
time, (T) 
frequency, 
RESULTS 


Shock Response Spectra.—As mentioned previously, response spectrum 
frequency the system. The response can represented the maximum 
value acceleration, velocity, displacement. example the shock re- 
sponse spectra used this study shown Fig. Frequency, plotted 
the abscissa; pseudovelocity, the ordinate; relative displacement, 
the scale sloping upward the left; and acceleration, the scale slop- 
ing upward the right. 

For undamped, linearly elastic, single-degree-of-freedom system in- 
itially rest and subjected base motion, the spectrum values response 
(that is, the absolute value the maximum quantity) can simply related (16). 
These relations also apply the spectrum values for the separate modes 
multi-degree-of-freedom system. However, when the total response the 
system calculated bythe simultaneous solution the governing differential 
equations motion, these relations are longer applicable. this study, 
the maximum value the total relative displacement was computed for par- 
ticular spring under certain conditions damping and applied forces, and the 
absolute magnitude was plottedas the spectrum value relative displacement, 
Therefore, this value the only one with direct physical significance. 
The spectrum values pseudovelocity and acceleration are computed from 
the same formal relations applicable single-degree-of-freedom system. 
The terms “v” and “a” are called “spectrum value pseudovelocity” and 
“spectrum value pseudoacceleration,” respectively. 

Amplification the objective this study 
the response free-base system subjected force applied the base 
mass with the response the same system considered sys- 
tem subjected base motion assumed represent the effects the actual 
disturbing force, helpful compare these responses terms ratio 


a 
4 
4 
is 
- 


called the amplification factor, AF. This amplification factor defined 
the ratio the “erroneous” spectrum value displacement the fixed-base 
system the “true” spectrum value displacement for the free-base sys- 
tem. separate amplification factor computed for each spring and 
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given subscript indicate the particular spring question. Thus, 
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General Form Results.—Thirty-four different problem groups were 
analyzed order study the effects the various parameters. Within each 
problem group, average twenty spring-mass systems were analyzed 
both fixed-base systems and free-base systems. Thus, the response was cal- 
culated for over 1,300 separate problems. Both absolute damping 

were included, and the amount was varied for each series prob- 
lems analyzed particular group. Response spectra were plotted for each 
spring each system for every set assumed conditions. 
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FIG. FORM RESPONSE SPECTRA 


Figs. and show the base displacement record that was obtained one 
the free-base systems for subsequent use the base motion input the 
fixed-base systems. These records are typical the others used for this 
particular forcing function. 

Figs. and represent typical response spectra that were obtained. 
Fig. the spectra for the first spring the systems with the amount 
damping indicated and subjected applied force, Fig. the spectra 
for the second spring the same systems under the same conditions. 
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indicates the general form the response spectra. The solid circles 
indicate the response the various free-base systems tothe particular force 
function under study. The open circles indicate the response the various 
fixed-base systems base motion that was recorded during the application 
the applied force particular spring-mass system. The response 
each system plotted function the first fixed-base natural frequency 
the system, The value for the free-base system which the record 
the base motion was obtained indicated onthe frequency scale. The free- 
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base frequencies, fj, each particular system plotted are interest and 
subsequent spectra these four values are also frequency scale. 

The free-base response and the fixed-base response are identical, except 
coincidentally, only for the particular system the base motion record 
was obtained. For all systems with natural frequency higher than this sys- 
tem, the fixed-base response larger than the “true” free-base response. 
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For all systems with natural frequency lower than this system, there 
random distribution values, although general the fixed-base response 
smaller than the “true” free-base response. The value the amplification 
factor varies depending the parameters specific problem; however, 
usually the fixed-base response the higher frequency range, that is, fre- 
quencies higher than the frequency the system which the base motion 
was recorded, much more error than the fixed-base response the 
lower frequency range. 
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The was computed forall the systems studiedand some these values 
are shown function frequency Figs. through 15. The difference 
the value the inthe lower frequency range and inthe higher frequency 
range can readily observed from these figures. The tends reacha 
peak value each the free-base frequencies the system which the 
base motion was recorded. However, the important observation that can 
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made from these figures that the remains large throughout the entire 
range higher frequencies. 

summary all the amplication factors computed for the first spring, 
AFj, presented Table The range values and the mean value are 
listed separately for the higher frequency range and the lower frequency range. 
Because there appears transition range between the first fixed-base 
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frequency the system which the base motion was recorded and the first 
free-base frequency that system, the values this range have not 
been included when computing the averages. The lower frequency range in- 
cludes all frequencies smaller than the first fixed-base frequency and the 
higher frequency range includes the first free-base frequency and all larger 
values. The mean value may slightly biased the higher range, because 
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more systems were studied the area which the amplification factor ap- 


ANALYSIS RESULTS 


The response the magnitude the effects the various para- 
meters; however, the amplification factor, that the ratio the responses, 
gives qualitative representation the effects. Whereas only values 

are presented this paper, indicate that were observed 
the other amplification factors. 

Effect Damping.—Viscous damping proportional absolute velocity and 
proportional relative velocity was studied. Figs. and indicate the 
typical results obtained. These figures are for the systems subjected 
applied force function 

Fig. 11, which shows the for absolute damping, indicates that increased 
damping substantially reduces the the region the first free-base fre- 
quency, This would expected because the damping was taken pro- 
portional the first natural frequency. the frequency increases above 
less. the region frequencies below fixed-base frequency, 
there only slight difference AF. Except for the region around 
the general shape the amplification curves similar. 

For the conditions shown Fig. 11, the average values for the higher 
frequency range were: 3.24, 2.37; 3.16, 2.31; AF3 3.21, 2.26; 
and 3.15, 2.24, where the first value listed for the smaller amount 
damping and the second for the larger amount damping. the lower 
frequency range the values were: 0.89, 0.79; 0.93, 0.84; AF3 
1.03, 0.83; and 1.01, 0.81. This ratio difference values gen- 
erally true the other problems which the systems and the applied force 
function were varied. Thus, doubling the amount damping the base mass 
and increasing five the amount damping the other masses decreases 
the about 40% the higher frequency range and about 15% the 
lower frequency range. The only major deviation from this pattern was found 
the systems which have large base mass. the systems the ap- 
pears insensitive increased amounts damping. 

Fig. 12, which shows the for relative damping, indicates that increased 
amounts relative damping have less effect the than did increased 
amounts absolute damping. The larger amount relative damping de- 
creases the about 15% the higher frequency range and 10% 
the lower frequency range. Again, the for the systems appears 
insensitive the amount damping. 

comparison can also made between the for absolute damping 
system versus the for relative damping the same system. Using the 
conditions given Figs. and and comparing values for the smaller 
higher frequency range 4.21; 3.16, 4.15; AF3 3.21, 4.21; 
and AF4 3.15, 4.05, where the first value for absolute damping and the 
second for relative damping. Thus, the same percentage damping 
used, relative damping will give slightly higher values than absolute 
damping would give. 
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Effect Applied Force Function.—A comparison for the various 
force functions applied the systems with relative damping 
0.02) shown Fig. 13. All the amplification curves show the same gen- 
eral tendency peak the four free-base natural frequencies the system 
which the base motion was recorded is, however, the val- 
ues these peaks and the shape the curves between peaks vary consider- 
ably. Considering the higher frequency range (that is, and listing the 
for Pg, P3, and that order, the average values are: 4.21, 
2.91, 3.57, 4.37; 4.15, 2.94, 3.28, 4.24; 4.21, 2.54, 2.98, 3.95; 
and 4.05, 2.33, 2.70, 3.76. The values the lower frequency range 
(that are more closely related and follow approximately the same 
distribution regardless the applied force function. 

Applied force reaches maximum amplification effect free- 
base frequency, and then produces almost constant over the en- 
tire range higher frequencies. Applied force reaches maximum effect 
the second free-base frequency, while produces large ampli- 
fication both the first and fourth free-base frequencies and 
Applied force produces large effects both and second free-base 
frequencies and Also, the peaks are much more pronounced 
for applied forces and 

Effect Magnitude Base Mass.—The effect the magnitude the base 
mass presented Fig. which shows the amplification curves for force 
applied the and systems with relative damping 0.20, 
=0.10). The systems have abase mass equal tothe sum ofthe story masses, 
the systems have base mass equal one the story masses, and the 
systems have base mass equal ten times one the story masses. 
would expected, the amplification several times greater for the systems 
with small base mass than for the systems with large base mass. 

The amplification curve for the system shows more ofatendency peak 
frequencies the system which the base displacement 
was recorded and than the other curves. the higher fre- 
quency range, considering AF3, and AF4 compositely, the for 
the systems varies from 13.50 5.10 with average about 7.6. The 
values for the and systems thehigher frequency range are more 
uniformly distributed with average about 3.5 and 2.2, respectively. 
Again, the values were more closely related the lower frequency range 
with composite average for the and systems approximately 1.0, 
0.8, and 0.8, respectively. 

This limited study indicates that approximate relation exists between 
the average amplification expected the higher frequency range and the 
magnitude the base mass. This relation indicates that the amplification 
factor varies roughly the square root the inverse ratio the magnitude 
the base masses. Thus, 


This relation crude because the amplification far from uniform, but 


does serve give some insight into the effect the base mass upon the er- 
rors that can expected. similar relation derived relating the maxi- 


on 
e 
aes 
‘= 
E 
| 
“he 
a 


October, 1961 


APPLIED 


0.10, 0.02 


FORCES, DAMPING 


© 
: 
=, 
a 
w 
5 On 
Ww 
e 
- 
° 


INTERACTION 


FREQUENCY, 


0.10 


0.20, 


FIG. 14.-AMPLIFICATION FACTOR FOR SYSTEMS, 
APPLIED FORCE DAMPING 


° 
Ab 


SYSTEMS, 
0.10, 0.02 


FIG, 15.-AMPLIFICATION FACTOR FOR 
APPLIED FORCE DAMPING 


O 


INTERACTION 
mum values obtained this study follows: 


Effect Spring Stiffness Distribution.—In all the problems 
previously examined, the spring-mass systems studied had uniform distri- 
bution mass and spring stiffness. Fig. shows the amplification curves 
for and systems with relative damping (Bo 0.02) subjected 
applied force Py. The systems have equal masses and equal spring stiff- 
nesses whereas the systems have the upper two-story masses and spring 
stiffnesses equal one-half the lower story values. The tendency peak 
the free-base natural frequencies the system which the base motion 
was recorded again evident and more pronounced the systems. The 
average values the higher frequency range for the and systems 
3.21, respectively. The amplification factor the lower frequency range 
varied approximately the same other problems that have been 
described. 

The base mass both the and systems equaltothe sum the story 
masses, but because the upper two-story masses are smaller the sys- 
tems, the base mass the systems smaller than the base mass the 
systems. Therefore, the free-base natural frequencies the systems are 
larger than those the comparison the amplification curve 
the systems subjected applied force (Fig. 15) and the amplification 
curve the systems subjected applied force (Fig. 13) shows re- 
markable similarity shape and tendency reach maximum value 
the second free-base natural frequency and Thus, was men- 
tioned before, there appears tobe evidence ofarelation between the frequency 
the applied force function and the free-base frequencies the system. 

Applying the relation derived previously for the effect the magnitude 
the base mass the average amplification the higher frequency range, 
was found that this case the relation was not valid. matter fact, 
the square root the ratio base masses moreapplicable than the square 
root the inverse ratio. explanation can presently given for this dis- 
crepancy. However, the maximum values are used, was found that 
the 0.7 power the inverse base-mass ratio still applicable. 

interesting compare these values the higher frequency 
range for the moderate-sized base-mass systems and Dwith the values re- 
ported Tung and Newmark (46). From the analysis ten-story building, 
they found values story shears from about 1.2 times larger than the 
design shears computed the procedure recommended the Joint Commit- 
tee Lateral Forces (1). Depending what system the base motion record 
was obtained for subsequent use input disturbance, the values for 
the and systems averaged from about only coincidence that 
these values fall the same range; however, the effect structure and foun- 
dation interaction may explain the reason for part the discrepancy between 
the results analyses and the behavior buildings subjected earthquakes. 


The purpose this investigation was study the errors that may occur 
when the foundation motion particular structure used determine the 
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response other dissimilar structures. common practice dynamic 
design problems use response spectrum concept todescribe specify 
the shock that the designed structure must withstand; however, developing 
this spectrum the base motion structure that has been subjected toa 
similar shock used assuming that this base motion accurate de- 
scription the shock. That say, the assumption made that the motion 
structure does not affect the foundation motion. the errors this as- 
sumption are large, then design based this concept will inaccurate. 
The designed structure may greatly overdesigned may quite in- 
capable withstanding the anticipated shock. 

this study various multi-degree-of-freedom, spring-mass models were 
developed represent dynamically actual structures. For all the model 
systems studied, was found that the lowest fixed-base natural frequency 
the model which the base motion record was obtained was important 
parameter. When the base motion recorded this model was applied 
base disturbance model with higher natural frequency, the 
response the fixed-base model always exceeded the “true” free-base re- 
sponse that model. When the base motion was applied models witha 
lower natural frequency, the fixed-base response, although variable, was gen- 
erally about equal the “true” free-base response that model. 

The magnitude the responses was found vary widely, depending the 
model parameters used. the higher frequency range, the fixed-base re- 
sponse was much times the free-base response. the lower fre- 
quency range, the fixed-base responses for the systems studied varied from 
maximum 2.59 minimum 0.25 times the free-base responses. 

Although there was large range the amount variance, important 
note that higher frequency systems always had larger fixed-base response 
and lower frequency systems usually had smaller fixed-base response when 
compared the true free-base response. Because the free-base and fixed- 
base responses only coincide for the exact system which the base motion 
was recorded, evident that errors will developed actual use. The 
actual dynamic properties and the behavior structure subjected shock 
can only estimated. Thus, when base motion record actual struc- 
ture used input disturbance for model the same structure, the 
calculated response will error unless the depicts the real 
structure. the model has higher frequency than the actual structure, the 
model analysis will predict responses that are larger than the actual re- 
sponses. the model has lower frequency than the real structure, the model 
analysis may predict smaller responses. 

Considering the model parameters used this study, the difference 
actualand computed response much less inthe range offrequencies smaller 
than the natural fixed-base frequency the system which the base disturb- 
ance was measured. Thus, would seem that base motion records obtained 
relatively high frequency structures, such stiff, short buildings 
equipment with high fundamental frequency, would provide more reliable 
data. Although this study primarily concerns buildings subjected ground 
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motion, the generalization the parameters makes applicable other 
shock transmission problems. 
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APPENDIX I.—NOTATIONS 


The symbols used the paper are presented here for ease reference 
and for the use discussers: 


alphabetic designation ofa series MDF spring-mass systems; 


amplification factor designating the ratio spectrum value 
displacement fixed-base system spectrum value dis- 
placement free-base system; 


spectrum value pseudoacceleration; 
alphabetic designation ofa series MDF spring-mass systems; 


alphabetic designation ofa series MDF spring-mass systems; 


coefficient viscous damping; 


coefficient critical damping; 


ie) 
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alphabetic designation ofa series MDF spring-mass systems; 
alphabetic designation ofa series MDF spring-mass systems; 
dimension force; 

undamped natural frequency (fixed-base system); 


undamped natural frequency (free-base system); 


time interval method; 


Subscript ith mass; 
Subscript jth natural frequency; 


spring stiffness; 
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dimension length; 


multi-degree-of-freedom system; 


mass; 

applied force function; 

undamped natural circular frequency (fixed-base system); 
undamped natural circular frequency (free-base system); 
proportion critical viscous damping; 

dimension time; 


shortest natural period vibration; 


time; 

relative displacement mass; 
relative velocity mass; 
spectrum value pseudovelocity; 


absolute displacement mass; 


absolute velocity mass; 


absolute acceleration mass; 


spectrum value relative displacement; 


parameter numerical integration method; 


critical viscous damping applied base mass 
absolute velocity); 


proportion critical viscous damping applied story mass 
(proportional relative velocity); and 


proportion critical viscous damping applied story mass 
(proportional absolute velocity). 
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ELASTIC RAYLEIGH WAVE EFFECTS DUE NUCLEAR BLASTS 


Melvin Baron,! ASCE and Charles 


SYNOPSIS 


nuclear burst the surface semi-infinite elastic half-space pro- 
duces atime decaying pressure pulse that acts over circular surface area 
increasing radius. This paper presents results for the stresses and dis- 
placements due Rayleigh waves that are produced such transient nor- 
mal surface pressure distribution. 

Expressions for the Rayleigh wave displacements and stresses produced 
concentrated load the surface the medium are used influence func- 
tions construct, means suitable integrations space and time, the 
corresponding quantities produced the transient normal pressure distribu- 
tions from nuclear numerical example presented which the 
Rayleigh stresses various depths and locations are given, The results were 
obtained using electronic computer, andthe numerical procedure de- 

discussion the numerical results and the applicability the theory 
actual situations presented, noted that the Rayleigh stresses pre- 
dicted the elastic theory are upper bound the stresses real me- 


Note.—Discussion open until March 1962. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the Engineering Mechanics Division, Proceedings the American Society Civil Engi- 
neers, Vol. 87, No. October, 1961. 

Assoc., Paul Weidlinger, Consulting Engineer, New York; and Adjunct Assoc. Prof. 
Civ. Engr., Columbia Univ. 

Chf., Programming Div., Army Ordnance Industrial Data Agency, Snowville, 
Va.; leave absence from the Mitre Corp. 
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dium. Their use for practical design purposes, such the design under- 
ground facilities, outlined. 


(Fig. 1), the major effects near the surface are produced Rayleigh Surface 
Waves. (These results have theoretically Rayleigh and are well 
confirmed seismological measurements.) This paper presents results for 
the stresses and displacements due Rayleigh Waves which are produced 
transient normal pressure distributions the surface the medium. Spe- 
cifically, time decaying pressure pulse acting over circular area in- 
creasing radius considered. 

The geometry the problem shown Fig. which the coordinates 
and locate the distance and depth the field point which the stresses, 
displacements, both, are obtained and the radius the disc 
pressure any time Expressions for the displacements and stresses due 


SURFACE 


FIG, 


Rayleigh Waves produced concentrated load with step pressure dis- 
tribution time, applied the point the surface the elastic half-spac 
have been obtained Chao, Bleich, ASCE, and Sackman. 
The corresponding stresses and displacements produced the transient nor- 
mal pressure which varies both space and time are constructed suitable 
integrations space and time the corresponding concentrated load quanti- 


ties. The results are valid points the medium for which the ratios and 
are small, that for points located shallow depth and sufficiently 


large distance from the radius the applied pressure disc. 

Numerical results are presented for normalpressure loadings produced 
thermonuclear explosions; specificallythe Rayleigh effects froma sur- 
face burst are evaluated. 


“Surface Waves Elastic C.C. Chao, H.H. Bleich and J.L. Sack- 
man, Report P-2066, Rand Corp., August, 
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Notation.— The letter symbols adopted for this paper are defined where 
they first appear, the illustrations the text, and are arranged alpha- 
betically, for convenience reference, the Appendix. 


MATHEMATICAL FORMULATION 


Expressions forthe stresses anddisplacements dueto Rayleigh Waves pro- 
duced transient normal pressure distributions the surface homo- 
geneous and isotropic linearly elastic half-space are derived this section, 
These results are obtained space and time the correspond- 
ing stresses and displacements due Rayleigh Wave effects from concen- 
trated, suddenly applied force the surface the half-space, 

The complete solution for stresses and displacements produced sud- 
denly-applied witha step distribution time has beentreated 
Pekeris and Lifson.4,5 The portion these quantities which are 
due effects were extracted fromthese results Chao, Bleich and 
the contribution certain poles the inversion integrals the 
Pekeris solution (the branch integrals appearing the complete solution can 
ignored). (The fact that surface effects can obtained the transform 
approach the contribution the complete solution was noted 
and Pekeris. This procedure based reasoning due 
Sommerfeld that the vicinity the Rayleigh wave front, the pole the 
leading term, exceeding all other contributions magnitude. 

The expressions for the Rayleigh stresses points located the coordi- 
nates (Fig. (1)) (measured from the instant surface pres- 
sure application the origin O), due concentrated suddenly applied load 
solution givenfor elastic material inwhich the equal, 


This relation holds for medium with Poisson’s ratio are3 


which (The notation and Stand for the real and imaginary portions 
the function 


Seismic Surface Pulse,” C.L. Pekeris, Proceedings, Natl. Academy 
Science, 1955, Vol. 41.7, 469. 

“Motion the Surface Uniformly Elastic Half-Space Produced Buried 
Pulse,” C.L. Pekeris and Lifson, Journal the Acoustical Society America, 
November, 1957, Vol. 29, No. 11, 

“Uber die Ausbreitung der Wellen der drahtlosen Telegraphie,” Sommer- 
feld, Annalen der Physik Folge, Vol. 28, 1909, 667. 
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denotes the shear modulus and refers the mass density the elastic 
medium, 


FIG, 


first step the evaluation the Rayleigh stresses due normal sur- 
face pressure distribution (r,t) the evaluation the corresponding stres- 
ses thatare produced bya surface pressure which constant magnitude, 


and uniformly distributed over circle radius From Fig. the 


stresses are obtained from the integral 


r+R 


(In the stress coefficient used Eq. andinthe subse- 


quent equations, the quantity the brackets gives the real part arguments 
the quantities and (Eqs. and 2f) which appear the equations for 


* 
(r,z,t).) Assuming that that is, that the results will valid 


points the medium which are located sufficiently large distances from the 
radius the pressure disc, the relation 


substituted into Eq. and, upon integration, 


The stresses produced normaltime-varying surface pressure dis- 
tribution (t), distributed over the circular disc radius obtained from 
the Duhamel integral 

(r,z,t) = p(t-7) aT. 608665 (6) 


which the Rayleigh stresses produced normal pressure with step 
function distribution time, obtained from Eq. Substituting Eq. into 
Eq. and simplifying, 
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Finally, the pressure p(t) expand radially time according 
the relation R(t). The Rayleigh stresses are obtained substituting 
R(t) into Eq. 


t R(7) 
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ec t-c T-yr-y 
s s 
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should noted that the expressions for the influence coefficients 
Eqs. are not valid the surface and hence, the integration procedure 
The determination the responses from transient distributed loads for 
requires separate approach starting fromthe transform expressions for the 
loads, equivalent procedure going the limit The reasons for 
this are explained the section for 
which are valid the surface are not presented this paper. 

The radial and longitudinal displacement components and uz, produced 
Rayleigh effects from the time decaying pressure acting over circular 
area increasing radius are obtained similar manner. The displace- 
ments u,* due concentrated surface load with step distribution 
time are, 


(8) 
2 Y 1 Re Z, Zo ( ) 
mur 


and 


which and are defined Eqs. and 2f. Proceeding through the 
successive space and time integrations(as the case for the stresses), the 
corresponding Rayleigh displacements and due the normal pressure 

p(t), which uniformly distributed over circular disc increasing radius 
Rit), are given the relations 


and 
R(7) 


(16) 


c 
* 
N 


and 


The evaluation the double integrals forthe stresses, Eq. 12, and the dis- 
placements, Eqs. 16and17, for given set pressure and range inputs, p(t) 
and R(t), using finite difference 
for both the space and time integrations. 


NUMERICAL PROCEDURE 


The evaluation the Rayleigh stresses and displacements particular 
field point (r,z) and time, from Eqs. 12, and respectively can per- 
formed numerically for given set pressure-time, p(t) and range-time, 
R(t), input functions. The numerical integration ofthe double integrals inthese 
equaticns represents major computational effort for cases practical in- 
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terest such normal surface pressure produced nuclear explosion. 
The variable upper and lower limits integration the variable add 
materially the numerical complexity problem. The numerical work re- 
quired the use large capacity electronic computer; 704 witha 
32K word memory was used obtain the results which will presented the 
following section. 

The present section describes some detail the numerical procedures which 
were used the analyses. The evaluation the Rayleigh stresses will 
illustrated. For convenience, the stresses are written the form: 


-9c 
which 
R(7) 


and and are defined as: 


and 


The integrals Eqs. and are integrated successively using Simpson’s 
Rule.’ The integration interval chosen for Eq. 19, thus yielding the 
pivotal values whichthe inner integrals areto evaluated. 

performed over even number subintervals. Each value 
multiplied the pressure form the integrand for the integration 
time Eq. 19. This integral turn evaluated Simpson’s Rule with 
47, thus yielding the stress the point defined the coordinates 
and the time The procedure repeated for each combination the input 
parameters r,z, andt. 

critical consideration from the viewpoint both accuracy and machine 
economy the determination optimum values for the integration intervals 
The computer program contains large number options which 
allow the printout key intermediate results such the inner integrals 


“Numerical Methods Engineering,” M.G. Salvadori and M.L. Baron, Prentice 
Hall, Inc., Englewood Cliffs, N.J., 1953. 
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thus permitting series preliminary sensitivity investigations 
made andoptimal values and beselected. These investigations 
are particularly necessary determining hg, due the nature the functions 
which appear the integral Figs. 3through show the 
the non-dimensional real part the argument and Zg. The extremely 
oscillatory character functions necessitate the ex- 
tremely fine grid spacing hg; some cases many 400 500 pivotal 
points were required. should also noted that the may neglected 
each case beyond value for the realpart the argument and 
Consequently, cutoff was incorporated into the program which set the 
equal zero for and greater than and less than -2; this re- 

sulted significant saving computation time. 

The input functions p(t) and R(t) which describe the transient normal sur- 
face pressure are chosen for typical nuclear explosions. Figs. and show 
curves for burst, used the computation the numerical results 
the following section. (The pressures produced the blast have been considered 
uniformly distributed overa circle radius R(t). This neglects small 
spike pressure, located near R(t), which early times contains only small 
the present study.8) each case, series sensitivity runs defineda cut- 
off time with corresponding pressure and range beyond which 
the pressure p(t) could taken zero negligible loss accuracy. 
The significance this cutoff that major Rayleigh effects are produced 
points the medium the relatively high intensity pressures the surface 
early times. 

The magnitude the numerical computations required that the integration 
Eq. made two steps. For the burst the integration was made 
used. 

Upon completion the sensitivity investigations which lead tothe determi- 
nation hj, the program was run without intermediate printouts. 
The computer program was prepared using the FORTRAN language. 


NUMERICAL EXAMPLE AND CONCLUSIONS 
APPLICATION NUCLEAR EXPLOSIONS 


infinite elastic medium bythe surface pressure from nuclear burst are 
evaluated the location 2000 ft, 3200 ft, and 5,000 from Ground Zero 
for various values the depth These values correspond the 10,000 
psi, 2000 psi, and 700 psi surface pressure contours for the explosion. 
The elastic constants for the medium were chosen follows: 


velocity shear waves 10,000 fps. 


“Cratering Megaton Surface Brode and R.L. Bjork, RM-2600, 
Rand Corp., June, “Theoretical Studies Ground Shock Phenomena, M.L. 
Baron, H.H. Bleich, and Weidlinger, Paper SR-19, Mitre Corp., October, 1960. 
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Poisson’s ratio 0.25 


mass density 5.186 


Figs. (a), (b), and (c) givetypical stress versus time histories the 
indicated locations for the Rayleigh stresses orr, and The time 
measured from the instant the explosion and hence, there 
delay before the first Rayleigh signals arrive the indicated field points. 
The results, which are typical all field points the medium, indicate that 
the stresses orr are larger than the other stresses that the importance 
the Rayleigh effects can generally judgedby Fig. shows 
the computedtime history forthe stress number depths the loca- 
tion 2000 ft. The rapid attenuation the Rayleigh stress with depth 
clearly shown. 

The peak values the stresses orr, and various depths can 
used obtain peakstress versus depth curves each location, defined the 
coordinate Figs. through show the rapid attenuation the peak stres- 
ses Orr, and with depth the location 2000 ft, 3200 and 5000 
respectively. each case, the peak values the shear stress were found 
comparableto the peak values stress Consequently, the 
stress versus depth curves for are not shown these figures. Attention 
again drawnto the fact thatthe formulas developed Chao, Bleichand Sack- 
and this paper cannot applied for points the surface, be- 
cause negative powers occur these expressions. Computational diffi- 
culties wouldarise for very the neighborhood and only 
results function the depth such that rather large stresses, and 
can expected near the surface (the vertical stress shear 
stress onthe free surface must vanish). The existence such high stres- 
ses inthe elastic analysis makes inapplicable actual situations near the 
surfacez and therefore the additional effort requires derive appro- 
priate analysis for shallow depth was not expended. 

interest investigate the effect change the shear wave ve- 
locity cg, that is, change the shear modulus the medium, the Ray- 
leigh stresses and Figs. and show the results which are 
6,000 fps but whose weight, 167 per cuft same asthe original 

medium. 

The peak stress versus depth curves 3200 ft, Figs. and can 
compared show the effect the decrease the shear modulus the 
medium. The peak stresses and each depth are considerably 

lower for the slower medium 6000 fps) than for the original medium 
(cg 10,000 fps). The characteristic rapid attenuation the peak stresses 
with depth shown for both media. 

Applicability the Results Actual Situations.—The analysis and the re- 
sults obtainedin this paper have been based onthe assumption that the medium 
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and remains linearly elastic under the application the surface pressures 
and that the effects from the slamming the bomb casing intothe medium are 
not considered. The elastic assumptions can violated two regions: (1) 
the crater region which surrounds Ground Zero, and (2) those points near 
the surface the medium which the Rayleigh stresses predicted elastic 
theory exceed the stresses that the rock materials under consideration are 
likely withstand elastically. 

the vicinity Ground Zero andthe surrounding crater region, the elastic 
theory not valid because, the extremely high pressures this region, the 
medium acts essentially fluid. The existence Rayleigh waves requires 
interaction between direct and shear stresses which does not occur ina 
fluid and one may conclude that the inception Rayleigh waves the crater 
region will inhibited. Because substantial portion the Rayleigh wave 
effects from the elastic theory originate from the high intensity pressures ac- 
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9.—CONTINUED 


ting the neighborhood Ground Zero, the results predicted the elastic 
theory willbe exaggeratedand the Rayleighstresses willtherefore upper 
bound the stresses that are produced the real medium. The problem 
the Rayleigh wave effects, which include the effect the crater, currently 
(1961) under investigation. 

The Rayleigh stresses predicted the present theory near the surface 
the elastic half-space willin some regions larger than the stresses that 
actual medium (such competent rock which could considered act 
elastically certain pressure levels) can withstand. Because the state 
stress surface waves inherently deviatory (that is, having considerable 
shear deformation), the medium would act inelastically such regions with 
consequent further reduction the Rayleigh stresses from the elastic theory. 
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(Dissipative effects will lead lower stresses due surface waves. Such 
waves viscoelastic medium have been considered Lubliner.9) 
Although the Rayleigh stresses predicted bythe elastic theory are anupper 
bound the true stresses medium, they are quite useful for practical 
design purposes. the results lies the applicationof the stress 
versus depth curves Figs. through and Fig. the design deep 
underground facilities. From these curves, possible determine each 
location those depths which the Rayleigh stresses will comparable to, 
less than, the stresses from the (dilatational) and (shear) waves which 
are also produced the surface pressures acting the elastic half-space. 
noted that each location 2000 ft, 3200 ft, and 5000 ft, the depths 
which Rayleigh effects become small compared and wave effects are 
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FIG. 13.—STRESSES DUE RAYLEIGH WAVES VERSUS DEPTH =5000 AND 
surface PSI 


relatively shallow due the very rapid attenuation the Rayleigh stresses 
with depth. This knowledge enables the designer use certain simplified 
obtain pressure inputs deep underground structures. such 
cases, diffraction information and shock spectra which are established the 
basis these simplified analyses could safely used for design purposes. 
detailed study these procedures the design underground structures 


“Surface Waves Viscoelastic Half-Space,” Lubliner, ONR Report Nr. 064- 
417, Report No. Columbia Univ., New York, April, 1960. 

“Theoretical Studies Ground Shock Phenomena,” M.L. Baron, H.H. Bleich, 
and Weidlinger, Paper No. SR-19, Mitre Corp., October, 1960, through 
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The displacements produced Rayleigh also considerable in- 
terest. Computations for these results are being run (as 1961) and are 


reported later time. 
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The symbols used the paper are for case reference and 
for the convenience discussers: 


y= 


velocity Rayleigh waves inthe elastic medium; 
velocity shear waves the elastic medium; 
integral See Eq. 20; 


time varying uniformly distributed pressure over 


concentrated load; 
range pressure distribution; see Fig. 


coordinates point the half-space; 
see Fig. 


coordinates; see Fig. 
time; 


displacements the and directions, respec- 
tively; 


influence coefficients for the displacements; see 
Eqs. and 18; 


specific weight elastic medium; 


functions; see Eqs. and 21b; 


ratio velocities; 


Poisson’s ratio; 
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shear modulus elastic medium; 
mass density elastic medium; 
stresses produced Rayleigh waves; 


influence coefficients for Rayleigh stresses; see 
Eqs. through 12; and 


direct and shear stresses produced Rayleigh 
waves. 
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INTERFACIAL MIXING STRATIFIED FLOW 


Enzo Oscar Macagno,! ASCE, and Hunter Rouse,2 ASCE 


SYNOPSIS 


effect state stratified flow inwhich the source interfacial 
turbulence would instability the interface itself, apparatus was devised 
which would permit the counterflow two superposed streams different 
densities maintained for considerable period time over appreci- 
able length uniform passage. means chronophotography injected 
droplets and analysis point samples fluid, measurements were made 
velocity and density distributions across the interface for various flow re- 
gimes. accordance with indications dimensional and stability analyses, 
results were presented the form dimensionless ratios interfacial shear 
and mass transport functions the Froude and Reynolds numbers. 


INTRODUCTION 


Density stratification occurs such natural fluid systems the oceans 
and the atmosphere well those that are man-made, and this can have 
important influence upon the flows that take place. many instances two 
more layers fluid different densities flow such way that they pre- 
serve their identity even whenthe fluids are miscible andappreciable diffusion 


Note.—Discussion open until March 1962. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Engineering Mechanics Division, Proceedings the 
American Society Civil Engineers, Vol. 87, No. October, 

Prof, Fluid Mechanics, State Univ. and Research Engr., Iowa Inst. 
Research, Iowa City, 

Prof, Fluid Mechanics, State Univ. Iowa, and Dir., Iowa Re- 
search, Iowa City, Iowa. 
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occurs across the transition zones between them. first approximation 
such layers can considered separated interface, rather than 
interlayer finite thickness, along which shear occurs and across which 
mass transfer takes place. 

The analysis the flow strata different densities requires knowledge 
conditions solid boundaries well along the internal boundaries be- 
tween the individual layers. The conditions solid boundaries are fairly well 
known, although the problem boundary roughness far from solved. The 
conditions inner boundaries, the contrary, have been only scantily ex- 
plored. sure, possible state them the case laminar flow, 
and some laminar-flow problems have been solved mathematically. However, 
soon the flow becomes unstable and the interface deformed waves, 
longer possible even the mean shear stresses. When tur- 
bulence develops, still another unknown introduced: the rate mixing 
the neighboring fluids and the resulting change the density distribution. 

many other problems, when mathematical analysis not yet feasible, 
one can resort studies model scale. However, model study can make 
ultimate sense without knowledge fundamental relationships which build 
either exact approximate criteria similitude. the present instance, 
approximate criteria forthe general case stratified flow must based upon 
the parameters derived forthe case discrete strata separated interlayers 
negligible thickness. 

The goal this investigation was experimental study interfacial in- 
stability and subsequent mixing due, however, conditions prevailing di- 
rectly the interface contrasted previous studies involving rather arbi- 
trary specialized boundary conditions. Only mention among those studies 
the ones that were the most useful during the present research, reference 
made the articles Taylor3 the stability superposed streams 
different densities, interfacial mixing and instability, 
and Raynaud6 underflow silt-laden water. All involved the flow 
one layer under over another that was stagnant except for the motion in- 
duced shear. 

Following the previously cited investigations, various papers have been pub- 
lished the 1950’s presenting new findings either the theoretical the 
experimental field. The series the hydrodynamical discussion stability 
was culminated the contributions Feldman7 and Graebel.8 Both 
papers contain solutions the Orr-Sommerfeld equation for two-layered strat- 
ified flows, taking into account the effect surface tension addition that 


“An Experiment the Stability Superposed Streams Fluid,” Taylor, 
Proceedings, Cambridge Phil. Soc., Vol, 23, 1927, pp. 730-731. 

“Interfacial Instability and Mixing Stratified Flow,” Keulegan, Journal 
Research NBS, Vol. 43, 1949, pp. 487-500, 

Harleman, Proceedings, Semicentennial Symposium Gravity Waves, NBS Circular 
521, 1952, pp. 

“Etudes des Courants Eau Boueuse dans les Retenues,” Raynaud, 4th 
Congress Large Dams, New Delhi, Question 14, 1951, pp. 137-161. 

the Hydrodynamic Stability Two Viscous Incompressible Fluids Parallel 
Uniform Shearing Motion,” Feldman, Journal Fluid Mechanics, Vol. Part 
1957, pp. 

“The Stability Stratified Flow,” Graebel, Journal Fluid Mechanics, 
Vol. Part 1960, pp. 
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density and viscosity. The case flow considered Graebel corresponds 
the one the present experiments, but his solution was confined essen- 
tially vertical flow, while that considered herein nearly horizontal. Grae- 
bel’s results indicate that rather low critical Reynolds numbers can ex- 
pected. the experimental side, the decade has seen frequent contributions 
Keulegan and his associates. Particular reference made the paper 
Lofquist,9 which presents results similar those for shear stress con- 
tained herein, although for dissimilar case flow. 
different and also somewhat abstract approach the problem inter- 
facial mixing was followed Rouse, ASCE, and experi- 
mental study the turbulent diffusion between miscible fluids different den- 
sities induced oscillation grid one the layers. The present in- 
vestigation, contradistinction, was intended cover process diffusion 
across interface due self-generated turbulence rather than turbulence 
imposed external means. Knowledge existing the time the present inves- 
tigation was started, 1957, indicated both the feasibility and the desirability 
this approach. The making the interface the plane maximum shear 
had already been under investigation the Iowa Institute Hydraulic Research 
Rouse and Yih. considerable amount work was still necessary 
reach the stage successfully controlling the flow, and develop technique 
observation which would reveal the mechanism interfacial mixing and al- 
low the researcher obtain quantitatively useful data. That the investigation 
was carried small scale was justified partly economical rea- 
sons but mainly previous experience showing that important relationships 
can often best revealed small scale which much more effective con- 
trol the phenomena canbe exercised. subsequent investigation should then 
continue the exploration the interfacial behavior much larger scale, ei- 
ther through controlled experiments the laboratory through the analysis 
field data obtained under natural conditions. 
Notation.—The letter symbols adopted for use inthis paper are defined where 
they first appear, the illustrations inthe text, and are arranged alphabeti- 
cally, for convenience reference, the Appendix. 


STABLE LAMINAR STRATIFIED FLOW 


The following analysis two-layered laminar stratified flow has the pur- 
pose defining the conditions under which nearly antisymmetrical flow 
possible. the case the uniform two-dimensional motion indicated Fig. 
the Navier-Stokes equations for the upper layer reduce 


and Stress Near Interface Between Stratified Liquids,” Lofquist, 
The Physics Fluids, 1960, pp. 158-175, 

“Diffusion Turbulente Travers une Discontinuité Densité,” Rouse and 
Dodu, Houille Blanche, 10, 1955, 522-532, 
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and for the lower layer 


The boundary conditions the top are 
and the bottom 


while the interface 


and 


Upon integration Eqs. and and use the conditions the interface, 
there results 


and 
where 
(15) 


The coefficients and are constants that can determined through the use 


pide 


and 

2 (41 hg + hy) ee 6 . 

Upon substitution and Eqs. and 11, follows that 


results 


and then follows that 


and 
(22 
the rates flow the two fluids are made equal, 
and 
which yields, finally, 
(25) 


‘ 

~ 


Evidently, the ratio kept nearly equal unity, the ratio will 
deviate still less from unity, and for all practical purposes the stratified flow 
can then considered consist two layers equal depth separated 
interface coincident with the plane zero velocity. 

Once possible produce the laminar flow two layers with velocity 
distribution which like sine curve, the stability the motion can stud- 
ied. With the inflection point the interface, there that disturbances can 
expected grow rather than any other level. this sense that an- 
tisymmetrical flow considered have inherent conditions instability, 
opposed other cases which disturbances develop elsewhere andare prop- 
agated the interfacial region. 


DISCUSSION INSTABILITY 


1945, concluded that, other things being equal, the relative sta- 
bility given flow should depend upon the velocity gradient, the wall distance, 


FIG, 1.—DEFINITION SKETCH FOR STRATIFIED FLOW WITHOUT MIXING 


and the kinematic viscosity, three variables which combine yield the nondi- 
mensional parameter 


(26) 


Evaluation this parameter for various points laminar flow near bound- 
ary, Shown Fig. results curve with definite maximum, the loca- 
tion which should correspond tothe zone minimum stability. Accordingly, 
comparison maximum values for different flows other things again being 
equal should permit their relative stability tobe assessed. Rouse also found 
that for conditions corresponding the lower critical Reynolds number for 
pipes (that which turbulence will longer persist) the maximum value 


General Stability Index for Flow near Plane Boundaries,” Rouse, Jour- 
nal the Aeronautical Sciences, 12, 1945, 329-332, 
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STRATIFIED FLOW 


was about 500 not only for pipes and wide channels but also for boundary 
layers and plane Couette flows. 

Evaluation the same parameter for the antisymmetric 
med Fig. indicates that two maxima exist: one the wall vicinity, be- 
fore, but another and much larger one the center line, where 

The latter, according this criterion, far the more susceptible dis- 
turbances comparable nature. If, moreover, the value 500 assumed 
torepresent the critical limit, this would found correspond Reynolds 
number only 83, which evidently full order below that for either 

wide channel the boundary region the antisymmetric case, for which 

2,250. Apart from the stabilizing effect gravity, which vanishes 

the density ratio approaches unity, the zone interfacial shear stratified 
fluid relatively susceptible instability the Reynolds type. 


FIG, PLANE POISEUILLE FLOW 


introduce the effect gravity behavior small distur- 
bances can examined through one-dimensional analysis12 interfacial 
waves. For each layer the stratified flow the following equations apply:13 


The two-dimensional case considered, because has been shown that there 
need take three-dimensional disturbances into account for the reason 
that they lead less rapidly instability (H. For long 


“Effect Turbulence and Channel Slope Translation Waves,” 
and Patterson, Journal Research NBS, 1953, 
“Advanced Mechanics Fluids,” Rouse, John Wiley and Sons, Inc., New 
York, 1959, 59. 
“On the Stability Three-Dimensional Disturbances Viscous Flow between 
Two-Dimensional Parallel Flows for Three-Dimensional Disturb- 


ances,” Yih, Quarterly Applied Mathematics, Vol, 12, No. 1955, pp, 434- 
435, 


waves the acceleration terms the y-direction can neglected. Additional 
will involve the neglect and and the 
hypothesis similar velocity profiles. After integration over across the 
strata, for each layer pair equations similar the following obtained: 


and 


Herein has been expressed terms the resistance coefficient, the 
ratio dy/U2h, and the pressure the interface. the linearization 


the equation performed and eliminated, one obtains the differential 
equation 


where the subscripts and correspond tothe upper and lower layer, respec- 
tively, and and indicate terms containing the resistance coefficients 
and 


and 


Eq. linearized form approximate equation for the interfacial 
waves. However, all the essential factors have been preserved it, that 
the pertinent dimensionless parameters involved can recognized: 


the Froude number; Ro, the Reynolds numbers each layer, included 
the corresponding resistance coefficients; velocity distribution co- 
efficients. should noted that two Reynolds numbers are equivalent one 
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STRATIFIED FLOW 


Reynolds number and the viscosity ratio. The ratio mean velocities 
not listed among the dimensionless parameters, because the imposed re- 
sional form, but new nondimensional variables will added. 


FIG, 3.—VARIATION ANTISYMMETRIC FLOW 


When the two layers are made equal thickness, the equation for the dis- 
turbance becomes simpler, and under the assumption equal resistance 
both layers with practically the same viscosity follows that 


which and are simpler expressions the terms and The re- 
lationship 


which satisfies this equation, represents disturbance which can made 
small desired assigning convenient values After substitution 
into Eq. equation for obtained, from which the condition that will 
make can derived that is, the condition for neutral disturbances 
which neither grow nor decay. Because this approximation the terms in- 
volving resistance laminar flow and for small differ- 
ences density, the rather simple condition 


. 
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Because the Froude and Reynolds numbers are related the case laminar 
flow the equation 


follows that the equation for the neutral curve the plane simply 


The question instability cannot, course, answered quantitatively 
two such partial approaches, even they are superposed some way. How- 


analysis 


FIG, 4.—SCHEMATIC REPRESENTATION THE STABILITY 
FUNCTION 


ever, one can expect acquire least qualitative indication from sucha 
combination. the lines represented the two analyses are plotted, zone 
the Froude-Reynolds plane determined (doubly hatched the Fig. 
which must certainly belong the stable region. least portion the 
neighboring zones inthe singly hatched regions shouldalso belong the stable 
portion the plane. all probability the boundary between the two zones 
determination the indicated function that the goal this investigation. 


EXPERIMENTAL APPROACH 


There simple and interesting experiment attributed Helmholtz 
which two fluids different density are made flow one over the other 
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carefully tilting the tube that contains them. The interface first seen 
very smooth, then regularly wavy, and finally disrupt with the on- 
set turbulent motion. The experiment canbe repeated indefinitely the two 
fluids are not miscible; can well with miscible fluids, 
although each repetition becomes more difficult. 


FIG, STAGES INTERFACIAL DISTURBANCE 


For the experiment with miscible fluids, care has tobe exercised tohave 
initially neat interface. regionclearly visible, dye canbe injected 
two fluids. The process mixing well two chem- 
ical indicators (such and sodium hydroxide) are used 
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such way that color will appear where mixing takes place. Figs. 5(a), (b), 
and (c) show interfacial conditions waves, with waves the verge 
breaking, and with turbulent mixing, respectively. The experiments were 
performed with fresh and salt water. Dye was used the first two, and chem- 
ical indicators were introduced the third. 

the Helmholtz experiment, which has been studied some its aspects 
Mittendorf,16 ASCE, when the difference densities small 
the symmetry the flow almost perfect and external disturbances can 
minimized; the phenomenon has transient character, however, and lasts only 
for minute two unless rather large scale used. Yet even large 
scale were adopted, every aspect the phenomenon would with 
time. Therefore, keep both layers motion and stabilize the occurrence 
different stages, apparatus was devised that would bring two streams 
flowing opposite directions into contact for reach considerable length. 


Observa tion 


FIG, 6.—SKETCH EXPERIMENTAL APPARATUS 


Fig. sketch givenof the apparatus used, which consisted mainly two 
tanks connected 36-in. conduit 2-in. rectangular cross section. 
third tank provided the necessary capacity for running the experiment dur- 
ing the time required establish the flow and perform the measurements. 
The apparatus was mounted tiltable frame order give adequate 
slope. Tank had lid that could closed during the experiments. 
Tank was open, and one its sides served overflow for the fresh wa- 
ter leaving the apparatus. 

perform experiment, the unit was first completely filled with fresh 
water and then salt water was started very slowly the bottom 
tank There, salt water was diverted radially and parallel the bottom 


“The Instability Stratified Flow,” Mittendorf, thesis presented the 
State Univ, Iowa, Iowa City, Iowa, 1961, partial fulfilment the requirements 
for the degree Master 
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plate and thence upward through wire screens. this way neat interface 
would formed, and the inflow could gradually increased the interface 
became more remote. When the interface reached the inlet the rectangular 
duct, the salt water would start flow down the slope toward tank and then 
Asthe salt water entered tank displaced fresh water, which 
moved through the apparatus inthe opposite direction and wasted over the weir 
tank The duration the experiments depended the rate inflow 
the salt water and the capacity tank 

During the experiment the discharge and the slope were varied sucha 
way that two layers equal thickness would form. aid this operation 
the two tailplates indicated Fig. had adjusted somewhat. One the 
main experimental difficulties was this short time and sucha way 
that appreciable disturbances could originate where the fresh and the salt 
water entered the working section and thus affect the flow within. very per- 
sistant disturbance another type was the intrusion the fluid interme- 
diate density resulting from mixing the two fluids; this, instead flowing 
away, would accumulate the end tanks, and sometimes had bled off 
order keep the interface clean. The great change the thickness ratio 
the two strata the tailplate sections tended stabilize the interface the 
approach section even while considerable turbulence was developing inthe zone 
under study. 


TECHNIQUES OBSERVATION AND MEASUREMENT 


addition measurement the discharge, slope, temperature, and geo- 
metrical dimensions involved, the velocity and density distributions the cen- 
ter line the central section the rectangular duct were subjected care- 
ful and detailed investigation. The choice this one section for the distribu- 
tion studies was determined for reasons not only convenience but achiev- 
ing state symmetry thetwo strata. was further assumed that the wall 
effects were small magnitude, that the central portion the flow could 
treated two-dimensional. Preliminary measurements well calcu- 
lations were made verify this assumption. 

The velocity measurement was based the technique photographing the 
motion small particles suspended the fluid, because the velocities were 
the low range and disturbances the flow were avoided. The well- 
known procedure chronophotography was adapted the case using 
timer Strobolume activated electrical clock combined with photocell. 
After considering and experimenting with several different tracers, the choice 
was made favor mixture two different fluids the proper proportion 
yield the desired density. For each experiment three mixtures n-butyl 
phthalate and xylene with bit white paint were prepared match the den- 
sities fresh water, salt water, and the intermediate mixture. The mixture 
was injected means hypodermic needles sucha way small 
clusters droplets essentially instantaneously but with practically distur- 
bance, was verified under conditions laminar flow. take pictures 
the droplets, narrow strip was illuminated with both continuous and inter- 
mittent light. The channel was made black plastic, sothat the droplets moved 
against dark background and left continuous traces the film with bright 
dots corresponding the flashes emitted the Strobolume. Usually some 
pictures were taken each run order ensure good average, es- 
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pecially the case turbulent motion. The traces left the film generally 
appeared have head and tail unequal length, and this helped deter- 
mine the direction the motion. effort was made synchronize the cam- 
era shutter and the flashing light sucha way make different heads and 
tails, because was found that mere chance was enough obtain the desired 
result most the pictures. 

The frequency the flashing light was controlled first means me- 
chanical contactor and then combination electric clock operating 
notched aluminum disk and electric circuit with photocell. beam 
light passing through the notches excited the photocell, which turn operated 
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the Strobolume. Both timers were built the Iowa Institute using standard 
equipment and instruments. 

Todetermine the density distribution, samples were taken different levels 
sufficiently small amounts minimize disturbance the flow and with- 
draw fluid only from the desired elevation. The samples were diluted after- 
ward larger volume and their electrical conductivity was measured de- 
termine the salinity and therefrom the density. Sampling and near the lat- 
eral walls was carried ina few experiments explore the conditions 
mixing that region. was found that the density distribution there was 
somewhat distorted with respect that the center and the results much less 
consistent. This vindicated earlier rejection optical methods based 
transverse beams light for determining the density distribution, because 


& 


they would include wall effects. was concluded that the use very fine hy- 
podermic needles would yield truer profile even they slightly disturbed the 
flow. 


DIRECTLY OBSERVED INTERFACIAL CONDITIONS 


Qualitative well quantitative observations can great value the 
study phenomenon which has modes occurrence yet disclosed. 
the naked eye and the photographic camera interesting features the strat- 
ified flow under consideration were made observable refraction due the 
density stratification through the addition dyes the injection parti- 
cles. The accompanying photographs illustrate the principal features 
thus observed. 

Fig. 7(a) the motion laminar with straight streamlines; the interface 
coincides with the line zerovelocity, which indicated droplets inter- 
mediate density appearing dots the center line. Fig. 7(b) corresponds 
interface with regular waves, and Fig. 7(c) interface with waves 
which break and start show irregularity and randomness. Fig. 7(d) depicts 
the motion after turbulence has developedand active mixing occurs across the 
interface. 

Based the qualitative conditions observed both visually and photograph- 
ically, each experiment was labeled with the symbol correspond- 
ing, respectively, uniform laminar motion, laminar motion with regular 
waves, irregular motion with breaking waves and incipient turbulence, and, 
finally, pronounced turbulence the interfacial region. These symbols were 
then plotted points inthe F-R plane. canbe seen Fig. four cor- 
responding regions are clearly distinguishable. express the correlation 
quantitative form, however, some measurable quantity had considered. 
The resistance coefficient, the shear stress the interface, and the rate 
mixing across were therefore evaluated described inthe following sections. 


THE RESISTANCE COEFFICIENT 


The stratified flow under investigation consists two streams flowing 
opposite directions and separated plane zero velocity. Considering the 
forces acting two fluid prisms length Ax, indicated Fig. one can 
write 


and 

where isthe wetted perimeter, the cross-sectionalarea eachstream, 
and (71)m and (72)m represent the corresponding mean shear stresses over 


floor, walls, and interface. From the two preceding equations, after elimina- 
tion one obtains 
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The masses and depend the density distribution and are tobe calcu- 


lated through integrals p(y) over the corresponding volumes. Approximate 
values and can obtained linear distribution density as- 
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sumed over the region where changes from the upper layer the 
lower layer. this way seen that 


and, similarly, that 


Substituting and the expression for yields 


where A/P the hydraulic radius. The resistance coefficient then 
given 


sin 


Since has treated the general problem, the solution for 
laminar flow ina rectangular pipe can evaluated therefrom. The corre- 


“Mémoire sur Influence des Frottements dans les Mouvements Réguliers des 


Fluides,” Journal Mathématiques Pures Appliquées, 2nd Series, 13, 1868, 
pp. 
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sponding resistance coefficient given 


Hence the ratio becomes 


This expression, although approximate one, shows that the correlation 
with the Reynolds and Froude numbers will essentially terms the 
relative thickness the mixing zone. 

The experimental determination the resistance coefficient required the 
measurement the discharge, the density distribution, the slope, and the di- 
mensions the cross the conduit. The corresponding values 
are given Fig. 10, which shows how they correlate with the Reynolds and 
Froude numbers each stream. There region the F-R plane where 
approaches unity, though inthe mean remains somewhat Where 
the flow isestablished stable, the value should exactly unity; 
this leads the conclusion that the lack establishment and the presence 
initial disturbances are responsible for being slightly greater than unity 
even the zone where stable laminar flow expected. That the mag- 
nitude the ratio otherwise steadily increases with and anticipated, 
quite evident from the diagram. 


DIMENSIONAL ANALYSIS FOR SHEAR STRESS AND RATE MIXING 


the interfacial shear stress and the rate transfer mass across 
the interface, these were next correlated with the independent variables the 
problem. The quantities and (Fig. 11) were selected represent the 
geometric and kinematic variables the interfacial layer, order divorce 
the analysis far possible from the dimensions the test section itself. 
Because the variations density and viscosity were small, the mean values 
and were considered accurate enough for the expression inertial and vis- 
cous forces. the other hand, gravitational effects being essentially depend- 
ent the difference between the two specific weights rather than upon their 
absolute values, was the corresponding independent variable. 
The dimensionless variables found through the approximate analysis anear- 
lier section are thus reduced Froude number and Reynolds number, the 
Instead considering directly the interfacial shear stress its ratio 
the viscous stress which also depends the above-listed variables, 
will analyzed. For turbulent flow the shear stress has been represented 
the sum two terms pu' Inthe present case values the 
shear stress greater than can expectedas soon the flow ceases 
correspond the steady-state solution the parallel laminar flow, though 
the initial increase not due turbulent transfer momentum much 
difference regime. desired, the ratio can considered “total” 
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viscosity coefficient. If, then 


dimensional analysis will immediately yield 
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For the rate transfer mass represented the mean velocity 
transverse motion across the interface, 


for which again follows dimensional analysis that 


Note should taken the assumptions made relative two different as- 
the analysis. One implicit inthe usual method selecting variables 
primary influence and disregarding those secondary The other 
involves the arbitrary separation regions flow whichare supposed 
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their own mechanical inner laws least first degree approximation. 
This certainly questionable step which should carefully checked under 
different Inthe present experiments was subjected verification 
changing the variable the ratio with satisfactory results. 


SHEAR STRESS THE INTERFACE 


The resistance coefficient gives only indirect measure the inter- 
facial conditions, because averages effects over the whole perimeter each 
stream, including the walls. thus variable which corresponds given 
rectangular section, whereas order obtain from these experiments pic- 
ture that focused the interface itself, desirable interpret them 
manner that minimizes boundary effects. 

The calculation the shear stress the interface for two-dimensional con- 
ditions similar that the resistance coefficient except that requires 
consideration the local value the shear stress and bottom bound- 
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STRATIFIED FLOW 


aries the conduit. involve considerable difficulty were not that 
the rather low values the Reynolds number ensured the existence laminar 
motion those boundaries (as could seen from the photographic records) 
even when turbulence occurred the interface. the equations for two ele- 
mentary prisms each written change momentum neg- 
lected, found that 


where the shear stress the interface. linear variation the density 
has been assumed over the thickness the mixing zone. After elimination 
the expression for becomes 


The flow the wall boundary layers being laminar, follows that 


y=-h 


accordance with the foregoing section, the ratio has been 
chosen represent the interfacial shear dimensionless form. Like the 
ratio must unity for parallel laminar should become greater 
than unity soon asthe flow becomes either nonuniform unsteady. order 
correlate the ratio with the Froude and Reynolds numbers, the latter pa- 
rameters were evaluated that only the region flow the vicinity inter- 
face need considered. The representative length chosen was the dis- 
tance betweenthe two maxima the velocity distribution curve, and the repre- 
sentative velocity was selected the maximum relative velocity be- 
tween the two layers. (This amounts, should remarked, treating the 
interfacial region layer with its own inner laws, approach which must 
tested further verify what extent may hold with different velocity 
distributions the rest the fluid.) 

The correlation with the variables 


and hence 
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shown Fig. 12. The scattering this case greater than the case 
the resistance coefficient, which probably due the difficulties involved 
determining exact values du/dy and Nevertheless, there certainly 
definite trend, and tentative curves have been fitted which show the approximate 
dependence the shear-stress ratio upon the special Froude and Reynolds 
numbers. These results are rather similar those reported Lofquist for 
another case flow.9 


TRANSFER MASS 


The onset instability brings about more active interfacial diffusion, 
which that point was due only molecular motion. The more turbulent 
the motion becomes, the more intense must the transfer mass from the 
denser the less dense fluid. The mean value the velocity mass transfer 
was calculated the following manner: the basis symmetry was as- 
sumed that, whereas across suchan interface the volumetric flux must zero, 
the mean velocity The longitudinal flux mass the entrance 
section each layer given simply the product the rate flow and the 
corresponding density. This flux increased with longitudinal distance for 
the fresh-water stream and decreased for the salt-water stream diffusion 
across the interface. The amounts increase and decrease can evaluated 
with respect the mass flux the central section through integrals the 
form u(y) over the corresponding section. approximate expres- 
sion these integrals obtained assuming that for the mixing zone 


which consists substituting for the density and velocity curves the tangents 
their inflection points. then follows that 
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Combination these two expressions yields 


and finally, 
as, = 12 L Ao 6.6.4.4. 8 6 (67) 


This velocity was correlated with the same variables was the shear 
stress, and through dimensional analysis was found that 


plot the plane 13) shows that the interfacial mixing 
varies systematically with these variables. Tentative curves fitted the ex- 
perimental points show trend similar that the curves for shear stress. 
This could certainly expected, since the shear stress depends diffusion 
momentum, while the rate measures the diffusion matter. 


CONCLUSIONS 


study has disclosed number significant facts regarding 
interfacial conditions stratified flow; these canbe summarizedas follows: 


Four successive regimes have been found exist the interfacial re- 
gion: 

The initial uniform laminar regime confined rather small zone 
the F-R plane, the latter parameters being defined terms variables 
describing layer the neighborhood the interface. 

Such uniform laminar motion first affected the onset motion 
non-turbulent type with relatively little secondary flow anda very small 
amount mixing. 

athird regionof the interfacial conditions are char- 
acterized breaking waves, which produce more effective mixing. 

Finally, the waves become complex spectrum with random char- 
acteristics, and turbulence develops the neighborhood the interface. 
The interfacial shear stress, expressed terms the ratio the ac- 

tual stress tothe viscous stress for the mean flow, shows adefinite correlation 
with the Froude and Reynolds numbers. 

similar correlation exists for the rate mixing across the interface. 
This anticipated, since thoughdetermined independently the shear 
measure momentum transport and the mixing measure masstrans- 
port. 
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general significance, these studies must now extended 
large-scale conditions the boundaries are remote from the interfacial 
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The following symbols have been adopted for use this paper: 


Froude number; 
resistance coefficient; 
thickness fluid stream; 
mass; 
wetted perimeter; 
pressure intensity; 
Reynolds number; 
mean velocity flow; 
differential velocity interlayer; 
mean velocity mass transfer; 
weight; 
interlayer thickness; 


one 
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dynamic viscosity; 

kinematic viscosity; 

mass density; 

normal stress; 

shear stress; 

mean shear stress along wetted perimeter; and 


stability parameter. 
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DYNAMIC RESPONSE BEAMS TRAVERSED TWO-AXLE 


Closure Robert Wen 


ROBERT WEN, writer wishes thank Enno Penno 
for his discussion. The greater influence higher frequencies shown the 
records the support reactions thetest bridge is, course, not unexpected/ 
The shearis the third derivative the deflection function while 
the bending moment the second. For the form the deflection function 
(Eq. 5), the contribution the shear and moment response of, say, the nth 
mode proportional 1/n and respectively. Furthermore, whereas the 
second mode does not contribute anything the bending moment mid-span, 
does affect the end reactions. 

Penno stated that the magnitude the total dynamic increment the re- 
actions could predicted computations within anarrow margin. From 
this might conjectured that the contributions the higher modes response 
would quite smallas comparedto that ofthe first mode, although their pres- 
ence was distinctly 

When highway vehicle moves over uneven bridge deck, essentially 
subjectedto external forcing function. this case, the importance damp- 
ing the response the vehicle and that the bridge would increase with 
the number energy input, assuming that reso- 
nant condition met. the number depending on, for example, the spacing 
the 3/4 in. in. boards used the tests appreciable, would seem that 
one should able observe reasonably distinct phenomenon resonance 
exhibited the bridge response. 

the other hand would seem more difficult observe field experi- 
ments the theoretically predicted “resonant speeds” due the timing the 
two axle loads applied (theoretically) smooth deck. this case essen- 
tially there are only two cycles energy input. This limits the magnitude 
the effects this “resonance.” Such limited effects could then over- 
ridden (and become undiscernible experiment) the effects other fac- 
tors that may have existed the actual test bridge-vehicle system, such 
the nonlinearity the vehicle springs and the deck unevenness, but have not 
been accounted for the theoretical calculations. 


October, 1960, Robert Wen Paper 2624). 
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BEARING CAPACITY FLOATING ICE 
Discussion Sebastyan and Penner 


and This paper was interesting be- 
cause the similarity the problem carrying capacity portland 
cement concrete pavements. 

The Canadian Department Transport’s responsibility design and 
evaluate the load carrying capacity airport pavements. For some time 
was felt that there definite need the conventional method 
design (Westergaard analysis) portland cement concrete pavements the 
basis field performance for those constructed inaccordance with Department 
Transport design standards and specifications, 

the summer 1959 and 1960, comprehensive load testing program was 
conducted number airport sites with varying subgrade soil condition, 
sub-base, base and portland cement concrete slab thickness. all, than 
fifty plate load tests were performed. 

One static load test was performed afree corner the slab withthe cor- 
ner freed from the adjacent slabs cutting through the slab with diamond 
saw. completion the test, the concrete was removed from area suf- 
ficient toclass the test and standard statictest was performed 
the top the base diameter plate was the minimum di- 
ameter that was used for the static base tests. One static load test was per- 
formed the diagonally opposite corner the removing any load 
transfer devices. additional static load test was performed the center 
the adjoining slab parallel thecenter line the runway. All thesetests 
were performed sound concrete slabs free cracks. Tested slabs were 
located minimum one panel distance from the edge the pavement. 

The deflection measurements for the corner tests were made means 
five dial micrometers placed line along the bisector the corner, and 
spaced approximately in. center. The first dial was placed close 
the corner possible with all deflections taken from independent datum. 

The deflection measurements for the center tests were made means 
four micrometers placed line along the bisector the opposite corners. 
Two dial micrometers were placed the plate, the other two in. from the 
center the plate. Absolute deflections were measured. 

load 2,000 was applied for sec obtain proper seating. The dial 
micrometers were thenzeroed. The first loadof kips was applied for min. 
Dial readings were takenat min, min, min, min, and min the two 
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dials closest tothe plate and min the three remaining dials. This pro- 
cedure was repeated with kip load increments until the pavement failed 
the capacity the load testing equipment was reached. 

note was made the load when the first crack appeared, and the actual 
load sudden loss bearing during load application. 

For free corner tests any possible edge support was removed for the fol- 
lowing distance from the corner: for in. concrete slab thickness 
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and for concrete slabs more than in. After completion the 
free corner test, three beams in, in. in. were sawed out. The 
flextural strengths (American Society for Testing Materials ASTM C78-59) 
the Young’s modulus the concrete, and the unit weight was determined to- 
gether with the compressive strengths (ASTM 
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After the free corner test was completed, all fractured concrete was re- 
moved and one static load test was performed the top the base course. 
Loading base test continued until 0.5 in. deflection was obtainedor 
until the capacity the test unit was After completing the static 
base test, three field California Bearing Ratio (CBR) tests and two field den- 
sity tests (ASTM D1556-58T) were performed. 


TIME MINUTES 


INCHES 


= 
z 
w 
w 
w 
a 
< 
= 


DEFLECTION RDG. 


FAILURE @) 


TIME VERSUS FIG, DEFLECTED SHAPE 
DEFLECTION CURVES SLAB UNDER SUCCESSIVE 
UNDER SUCCESSIVE LOAD LOAD INCREMENTS 
INCREMENTS 


Load testing rigid pavement slabs was conducted when the corners and 
edges the slab were curled downward. Testing commenced between a.m. 
and p.m. local standard time. Slabs unshaded until the 
time the test. Tests were completed after sunset. 


var 


standardauger hole was drilled adepth top the slab 
each free corner test location. Soil samples recovered were classified 
accordance with the Corps Engineers soil classification system, and 
the moisture content was determined intervals. Sketches and photo- 
graphs were made the crack pattern for the slabs tested destruction. 

Typical examples thetype data received from the field and the method 
presentation are given Figs. 13, 14, 15, 16, and 17. 

The load deflection and time deflection relationship, together with the 
shape the deflection bowland the modulus subgrade reaction were deter- 
mined various loadings and deformations. short summary the test re- 
sults given Table 

The tables give, addition the physical data, the comparison between 
the actual load carried and the load carrying capacity determined the basis 
the Westergaard equations and the Meyerhof analysis. 

The measured load carrying capacity was found considerably higher 
almost every instance under the test conditions used than that predicted 
the Westergaard equations. This shows that limit design method might well 
the correct approach. the other hand, the reason for the lower values 
obtained using the Westergaard equation might attributed the assumption 
made evaluating the influence the subgrade the load carrying capacity 
the system. 

Instead assuming the subgrade heavy liquid, the total subgrade 
reaction can determined the basis the measured deflection bowl, to- 
gether with the integration the effect the modulus subgrade reaction 
various deflections. 

considered that these test results are most significant. Together with 
the evaluation the effect traffic and alimit design method may 
form new basis for portland cement concrete pavement evaluation. The ob- 
ject this presentation make available field data which indicates the in- 
herent applicability Meyerhof’s portland cement concrete slabs. 

intended thatthe evaluationof these test results and the possible modi- 
fication the Canadian Department Transport’s rigid pavement evaluation 
method will the subject separate paper. 

Notes Table Load Test Results.—(The yield load, center loading con- 
dition.) some cases, the yield point could determined from break the 
load-deflection curve (Fig. 16). other cases, the load-deflection curve ex- 
hibited definite yield point (Fig. 17). these latter cases yield load 
reported the table, although some cases yielding must have occurred. 

Originally, the panels were 100 long, but the time the test, 
these were sub-divided one more transverse cracks. case, was the 
length the sound section less than the width. 

Anassumed “E” value was used. The “E” the mean oftests onthree 
beams from this airport. assumed flexural strength was used for these 
sites. 

The value reported the mean breaks from beams cut various lo- 
cations inthe field that concrete was approximately 
old. While the beams were not cut fromthe slabs actually tested, they were cut 
from similar concrete placed under the same conditions. Therefore, the average 
strengthof the beams should represent the average strengthof the 
The value reported has been computed from the strengthof 28-day construction 
control beams cast during the day when the test section was constructed. The 
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28-day strength was Portland Cement Association (PCA) 
specifications.24 The value reported was computed from the mean result 
repetitive plate bearing tests the subgrade, according the Department 
Transport Procedures. 

Procedures.—The value reported was computed from repetitive plate 
bearing test each location, according the Department Transport Pro- 
cedures. The value computedfrom the mean result four repeti- 
tive plate bearing the base this Airport, according tothe Department 
Transport Procedures. The value reported was determined from one test 
this airport. 

Keyed corners were considered free corners for the purposes these 
computations. realized that this conservative assumption. theo- 
reticalinvestigation the revealed there was some 
justification for assuming semi-infinite plate condition and this was done 
six cases airport This not conservative assumption, but was felt 
that was more realistic than disregarding the effect the dowell. 

The factor for use Fig. (a) was taken for each loading con- 


dition, which the measured “modulus rupture.” For corner loading 
(90° Wedge) equivalent was determined the radius quadrant whose 
area equaled the area the plate used for the test. For edge loading (semi- 
infinite equivalent was determined, the radius semi-circle 
whose area equaled the area the plate usedin the test. For each loading con- 
dition the equivalent was further modified the addition depth factor 


since all test locations <2h. 
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ANALYSIS NUMERICAL PROCEDURES* 


Closure Annabel Tong 


ANNABEL ASCE.—The writer wishes mention the pos- 
sibility programming the computation steps this numerical procedure 
highspeed digital computers. spite the simplifications used this meth- 
od, the computation work involved, done manually, still time-consuming 
especially for the bending buckling framed structures. The majority 
the work inthis analysis can programmed smallor medium sizedcom- 
puter with controlled degree accuracy. For example, the case column 
buckling, Eqs. 18, 19, 20, and are first derived for column specific 
material and cross-sectional shape. Corresponding values 12, and 
Table can computed the computer and stored inits memory unitsor di- 
rectly printed out table formats for future uses. The steps for constructing 
Figs. 13, 14, and can programmed and results can arranged groups 
under given loadings. For example, the various values andthe corresponding 
could pairs under given P/A andpre-selected deflec- 
tion The next step regroup these data storing and the corres- 
ponding pairs for given P/A and logical decision, the computer 
will then choose the maximum number among agroup numbers represented 
under each pair P/A and This step equivalent the step for de- 
termination maximum point curve shown Fig. 15. The chosen 
maximum length will furnish one the buckling curves plotted Fig. 
The program for this problem may written many separate steps. The 
results step wouldbe printed out andchecked before the next step 
accomplished. Once such program available, could used repeatedly 
for columns other cross-sectional shapes and materials with slight modi- 
fication. 

generalized column curve would possible for members various sec- 
tions and endconditions shape factors and effective length factors are intro- 
duced this analysis. basis for design, member rectangular section 
may assumed have shape factor columns hinged ends have 
effective length factor Then the factor for member section other 
than rectangle should appliedto the eccentricity ratio R2) and the fac- 
tor for restrainedcolumns should applied the slenderness ratio (L/R). 
The column are plotted scales P/A versus L/Rfor various 
The effective lengths restrained determined methods 
commonly used the elastic analysis. The determination the shape factor 
crude approximation. References for the values for various cross- 
sectional shapes are available (7). 
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RESPONSE MULTI-STORY STRUCTURES 
Discussion Richard Troy 


RICHARD TRoy, stimulating paper prompted the 
writer make some calculations for the four-story frame. These consisted 
estimating the first mode period the bare steel frame, ascertaining the 
approximate effect the usual solid concrete fireproofing columns and gir- 
ders the value the first mode period and comparing the seismic coef- 
ficients tabulated author Table with those obtained applying the 
provisions the latest edition the Los Angeles City Building 

estimating the period, the bare steel frame was loaded the floor levels 
with the loads shown Fig. lof the paper and the columnand girder moments 
were obtained the portal The story deflections were then computed 
and the total top story deflection was obtained. This value was substituted 
the formula 


which the total top story deflection inches and the coefficient 
the average the coefficients the formulas for the free period 
weighted cantilever loaded cantilever, the coefficients being 
0.32 and 0.26, respectively. value sec was compares 
with 2,01 sec obtained the author. 

The effect solidfireproofing was found calculating the increase stiff- 
ness the columns lightweight concrete with was assumed 
and the moments inertia the cracked and uncracked sections were found 
assuming rectangular concrete section with in. minimum cover over the 
steel sections, The uncracked fireproofing some cases increased the stiff- 
ness more than 100%. The increase, considering the fireproofing cracked, 
was less, was expected, assumed that the solid fireproofing 
increases the stiffness columns and girders 75%, then the total top story 
deflection will 57% that the bare steel frame and the free period will 
reduced 25%. This would givea period for the fireproof frame 1.92 
0.75 1.44 sec, that would equivalent tothe undamped structure with shear 
walls Berg’s paper. The witha smallern value would 
increase the stiffness the frame even more than light-weight concrete with 
resulting reduction free period. account was taken the contribution 
the floor the stiffness the girders with portion thereof acting 
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Any yielding slip the girder column connections, the other hand, 
would tend counteract this increase stiffness. 

The comparison seismic coefficients revealed some interesting differ- 
ences. 1960 the Los Angeles City Building Code was revised anda more 
elaborate set procedures was designed seismic loads. Previous 
this, the seismic coefficients for stories less were de- 


termined from the formula where was the number stories above 


4.5 
the one under consideration. The new code defines base shear 


which factor dependent the type building construction, isa 
factor 


dependent the period 


the structure and the total dead load the structure with qualifica- 
tion added for storage buildings. formula given for the period based 
the height and width direction motion the building. further 
stated that number stories for all buildings with moment re- 
sisting frame carrying 100%of the seismic loads, and that need not taken 
less than 0.10 

For the four-story frame the code would give 0.4 sec, value which 
the value given the author. Nevertheless, the use this value for 
will result seismic coefficients that are line with those given 1(a) 
the paper. The following the seismic coefficients according 
the 1960 Los Angeles City Code for the values listed: 


(Bare steel frame) (Fireproofed frame) (Los Angeles 
Code values) 


From this appears that one ill-advised substituting the actual free 
period flexible low building into the Los Angeles Code formula. also 
indicated that the procedures the code are device for obtaining seismic 
loadings that reflect past damage experience rather precise methods 
for finding the dynamic characteristics structures. 


re = 0.05 (6) 


DISCUSSION 


The seismic coefficients eight-story frame were also determined 
see they were line with Berg’s values. The results are follows: 


0.8 sec 
(Los Angeles 
Code value) 
0.065 
Ce = 0.056 
Cs = 0.052 
0.048 
Co = 0.040 
0.036 


evident that the new Los Angeles Code has reduced the seismic loadings 
appreciably over those the old code especially the upper con- 
sideration Berg’s data forthe viscous damped frames would seem justify 
this except for the Centro, Calif. ground motionas affects the four-story 
frame. would appear from his data also that higher seismic loadings than 
the code would prescribe should used for structures where little damping 
will present, even though their construction would justify the lower coef- 
ficients stated the code. 

With regardto damping non-structural elements, the writer feels 
that too much emphasis has been placed their furnishing such damping. The 
modern multi-story building does not often have tile masonry partitions. 
The usual practice use movable partitions that extend the ceiling 
line. Where plaster partitions are used, the plaster stopped the ceiling 
line unless fire rated wallis required, the partition plastered 
the underside the floor above. This the case around elevator shafts, 
stair wells and duct shafts. Where the plaster stopped the ceiling line 
the light gage metal studs extend beyond and are fastened the underside 
the floor above. These studs bending about their weak axis will furnish only 
small amount resistance motion. The fire rated walls therefore furnish 
the lion’s share non-structural damping together with the curtain wall, which, 
keeping with modern practice, lightweight and many cases largely 
result, the fire rated partitions will damaged earthquake they 
are expected furnish damping. There then arises the question whether such 
damage should tolerated some sort fire rated joint should devised 
the tops such partitions allowthe floor above move without dragging 
the partitions with it. The initial cost the fire rated partitions typical 
multi-story office building will considerable sum money. Need- 
less say the cost repairing replacing them they are badly cracked 
might well equal exceed their initial cost. Therefore, the decision count 
non-structural items for damping should reached after consideration 
given the economics involved. 

With regard damping provided inelastic deformation, the writer feels 
that Housner’s thesis might winacceptance for certaintypes structures such 


Pe 
é 
sip 
— 


October, 1961 


water towers and certain industrial doubtful many prac- 
ticing structural engineers building officials would accept for the design 
multi-story apartment office buildings. 

The author has made valuable contribution the literature seismic 
hoped that other studies can made actual buildings 
designed for seismic forces and that field measurements their dynamic 
characteristics during various stages construction available for use 
such 


102 
j 
A 


103 


COMMENTARY PLASTIC DESIGN 


Closure Committee 


CLOSURE WRC-ASCE COMMITTEEON PLASTICITY RELATED 
DESIGN.—This closure applies the seven progress reports which com- 
prise the Commentary Plastic Design Steel. The Joint Committee 
pleased note the several constructive discussions that were submitted. They 
serve indicate possible refinements the plastic method design and 
state areas where misconceptions may exist. 

The Joint Committee considers that has misinterpreted the contents 
and purpose the 

The discussion Art. 6.2 (Local Buckling) the Commentary 
Fialkow, ASCE, primarily concerned with the fact that Eq. 6.16 (which 
gives the maximum allowable beam-depth web-thickness ratio for members 
subjected axial bending moment) basedon maximum strain 
four times the yield strain. Fialkow shows with two examples that theoretically 
much larger strains can exist hinge thetime the failure mechanism has 

Fialkow’s critical review the fundamental limitations Eq. 6.16 ap- 
preciated. However, the members the Joint Committee responsible for the 
publication the Commentary feel that the local buckling guides set forth 
Art. are satisfactory. not statedin this article that the strain must 
limited Eq. 6.16 suggested guide for checking the critical d/w 
recommendation Ref. 6.27, (In ASCE Manual 41, This Ref, 6.25). 
can used.” The limitations Eq. 6.16 are also clearly stated Art. 6.2, 
and the guide recommended for structures for which ordinarily deflection 
and rotation analysis would performed. Furthermore, the curves Fig. 
6.11 enable the designer determine the maximum allowable value d/w 
when larger strains are encountered, such more complicated structures; 
also, such relationships should helpful providing bounds cases involving 
dynamic loading. 

The Joint Committee feels that Eq. 6.16 appropriate for plastic design 
for the following reasons: 


The development Fig. 6.11 which Eq. 6.16 basedcontains num- 
ber conservative assumptions (Ref. 6.27). Local buckling the web as- 
sumed take place soon the yield stress reached, the start the 
strain-hardening range. Thus, the material not assumed strain-harden, 


Published seven Progress Reports follows: Report and Chapters 


July 1959; Report No. and Chapters October 1959; Report No. Chapters 
Jan, 1960; Report No. and Chapters and April 1960. 
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however, strain-hardening actually does occur, especially when the moment 
not uniform. 

The theory onwhich the curves are based does not consider 
the large post-buckling strength which available even for thin 

shown 9.8 the Commentary that for more complex struc- 
tures major part the rotation the first hinge takes place after almost 
all the strength the structure has been realized (Fig. 9.9). Thus, even 
large strains are present the formationof the mechanism, the strains 95% 
the maximum load may considerably less. Therefore, undue concern 
about the exact magnitude not warranted, Also impossible 
obtain reliable estimate the strains failure from the hinge angles, 
these are determined under the assumption that the hinges form point. 
more exact strain calculation involves complicated elastic-plastic analysis; 
only for the simplest structures are such analyses expedient. 

quantity experimental evidence has shown that the lo- 
cal buckling rules 6.2 are satisfactory. Much this evidence given 
Ref, 6.27 whichexperiments were made substantiate the local buckling 
theory. Further evidence can seen Figs. 5.11 the Commentary 
where the full-scale frame test results are tabulated. may observed that 
the predictedultimate loadwas attained even fabricated from mem- 
bers which had high d/w ratios (for example d/w for the 10B17 members 
which were part severalof these frames). About forty experiments beam- 
columns have shown that the primary cause always lateral-torsional 
buckling general column instability. (The results several these tests 
are given Figs, 7.19, 7.20, and 7.21 the “Commentary Plastic Design 
Steel.”) Several these experiments were performed laterally braced 
8B13 columns (d/w 35). Even for 0.6 failure was due general 
column instability. 


The Plastic running account conducted research, 
and such cannot expectedto have the “last say” about many ofthe prob- 
lems studied it. The problem inelastic instability certainly compli- 
cated one, and hoped that future research toa better under- 
standing the problem. particular, the Commentary implies, future re- 
search should provide experimental data that can used for comparison with 
the theory which Eq. 6.16 was based. 

Hashim Hamzawi, ASCE, correctly stated that “critical” section 
might occur some intermediate point tapered haunch the adequacy 
the haunch checked only its two ends. However, Hamzawi did not say that 
the magnitude the inadequacy smallindeed, For the linear moment diagram 
given Hamzawi’s example can shown that the deficiency resisting 
plastic moment the critical section only from 1%to 5%of the total that 
(8.7) 

Using Hamzawi’s symbols, the following reduction moment occurs 


Aspects Concerning Inelastic Instability Steel Structures,” Thurlimann, 
Proceedings, ASCE, Vol. 86, No, ST1, January, 1960, 99. 
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For all practical purposes this amount can neglected design using 
adequate load factor recommended Art. 4.8 the Commentary. 

Garden has presented interesting refinement for including the effect 
spread plasticity the calculation deflections the slope-deflection 
method. The refinement seems more adapted orderly computation than the 
other precise methods the report. 

would interesting obtain values the commuted curvature 
gravity commuted curvature for average wide-flange and shapes. This 
would relatively simple triangular moment diagrams the type shown 
the discussion, the case members with constant applied moment equal 
along finite length and members with plastic hinge forming re- 
gion uniform load, additional limitations would have Pos- 
sibly consideration these limitations might involve introducing the effect 
strain-hardening into the calculations. 

The Joint Committee pleased note the rapid strides made the field 
structural design insteel from 1940-1960. large measure, these advances 
have resulted from understanding the relation between ductility and the 
performance steel structures. The Commentary suggests how the property 
ductility may used structural design and reviews the necessary secondary 
design considerations which insure that ductility fully and safely utilized, 
The resulting design method, known Plastic Design Steel, has been used 
many structures built the United States and abroad, 

The seven progress reports which constitute the Commentary have been as- 
sembled one volume June 1961as ASCE Manual Engi- 
neering Practice 41, addition the topics covered the progress re- 
ports, this extensive list references plus author 
ject indexes. 

The Joint Committee well aware that continuing researchon the inelastic 
behavior steel structures necessary. Some this work currently (1961) 
underway. Its completion will make possible utilize more fully the poten- 
tial plastic design such structures multi-story frames and stiffened 
plate assemblies, potential that experience has shown will combine effectively 
the dual objectives safety and economy steel structures. 
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